Effect of Long-lived Strongly Interacting Relic Particles on Big Bang Nucleosynthesis 
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It has been suggested that rehc long-lived strongly interacting massive particles (SIMPs, or X 
particles) existed in the early universe. We study effects of such long-lived unstable SIMPs on big 
bang nucleosynthesis (BBN) assuming that such particles existed during the BBN epoch, but then 
decayed long before they could be detected. The interaction strength between an X particle and a 
nucleon is assumed to be similar to that between nucleons. We then calculate BBN in the presence 
of the unstable neutral charged particles taking into account the capture of X'^ particles by 
nuclei to form X-nuclei. We also study the nuclear reactions and beta decays of X-nuclei. We find 
that SIMPs form bound states with normal nuclei during a relatively early epoch of BBN. This leads 
to the production of heavy elements which remain attached to them. Constraints on the abundance 
of particles during BBN are derived from observationally inferred limits on the primordial light 
element abundances. Particle models which predict long-lived colored particles with lifetimes longer 
than ~ 200 s are rejected based upon these constraints. 

PACS numbers: 26.35.-l-c, 95.35.-|-d, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

There has been considerable recent work on the elFects 
of decay or annihilation of exotic particles on light ele- 
ment abundances [l|, 0, 0, 0, H, B S Hij iS] • Since standard 
big bang nucleosynthesis (BBN) predictions of light el- 
ement abundances are more or less consistent with ob- 
servations, changes of abundances relative to those of 
the standard BBN cannot be large. This makes it pos- 
sible to constrain theories beyond the standard model 
through their consistency with observed light element 
abundances. Moreover the decay process of massive par- 
ticles might change the lithium abundances providing a 
solution to the lithium problems. Recent studies suggest 
that radiative decay could lead to the production of ^Li 
to the level at most ~ 10 times larger than that observed 
in metal-poor halo stars (MPHSs) when the decay life 
is of the order of ~ 10^ — 10^^ s [9], and the hadronic 
decay can be a solution of both the lithium problems al- 
though that case gives a somewhat elevated deuterium 
abundance [1, [1]. 

The possibility of the existence of heavy (to 1 GeV) 
long-lived color flavored particles has been discussed in 
scenarios of split supersymmetry [l^, |Tl[, and weak 
scale supersymmetry with a long-lived gluino [T^ or 
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squark |13| as the next-to-lightest supersymmetric par- 
ticles. Those heavy partons would be confined at tem- 
perature below the deconfinement temperature Tc ~ 
180 MeV inside exotic heavy hadrons, i.e., strongly in- 
teracting massive particles (SIMPs) [ij]. Under the as- 
sumption that the X particles are in statistical equilib- 
rium with the thermal background in the early universe, 
Kang et al. lj| estimated a relic abundance of those 



hadrons based upon a comparison between their anni- 
hilation rate and the Hubble expansion rate. In this way 
the estimated relic abundance can be written 
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where Nx is the number density of the X particles, and 
s = 27r^(7*sT'^/45 is the entropy density with ~ 10 the 
total number of effective massless degrees of freedom [lB| 
just below the QCD phase transition. T is the tempera- 
ture of the expanding universe, R is the effective radius 
for annihilation of the X particles (i? ~ GeV~*), Tb 
is the temperature at which the X-particles are formed, 
and TO is the mass (to ^ 1 GeV) of the heavy long- 
lived colored particles. This relic density corresponds to 
a number fraction of 
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where ni, is the number density of baryons. We therefore 
assume the existence of the long-lived heavy hadronic 
particle X and study effects of such particles on BBN. 

Experimental constraints on hypothetical SIMPs have 
been delineated in flU, [3] ■ The effect of new neu- 
tral stable hadrons on BBN was studied in [l^. They 
assumed that the strong force between a nucleon and a 
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stable hadron is similar to that between a nucleon and 
a A hyperon and that most new hadrons end up in a 
boimd state of *He plus the hadron after BBN. The re- 
sult of their analytical calculation showed that the stable 
hadrons would be preferentially locked into beryllium. 
In other words, beryllium has the largest fraction Ax /A 
of bound states with the hadrons among the light ele- 
ments produced in BBN, where the A and Ax represent 
a nuclide A and a bound state of A with a hadron X. 
Mohapatra and Teplitz [lO] estimated the cross section 
for an X to be captured by "^He and claimed that the 
fraction of hadronic X particles captured by ^He nuclei 
is smaller than that assumed in [l^. Therefore a large 
fraction of free X particles would not become bound into 
light nuclides. 

In this paper we carry out a consistent calculation of 
BBN in the presence of a hypothetical long-lived SIMP 
X^ of charge zero assuming that the X^ nucleon interac- 
tion is of a similar strength to that between two nucleons. 
Binding energies of bound states of X-particles and nu- 
clei, which we call X-nuclei, are estimated. Rates for X 
capture by nuclides, as well as the nuclear reactions and 
P decay rates of X-nuclei are also estimated. We cal- 
culate BBN including the particles as a new species 
taking account of many reactions related to X'^ particles 
in a network calculation, and study the effects of par- 
ticles on BBN. In Sec. HI] assumptions regarding the X 
particle, estimations for binding energies of X-nuclei and 
various reaction rates are described. In Sec. IIIII results of 
the network calculations are shown, and the constraints 
on parameters of the X*^ are derived from a compari- 
son with observed primordial light element abundances. 
Conclusions of this work are summarized in Sec IIVI 



II. MODEL 

We have added the X particles and relevant X-nuclei 
Ax as new species. Their reactions have been added 
to the BBN network code of Refs. [H [H. Nuclear re- 
action rates for the standard BBN have been replaced 
with new rates published in Ref. 12311 and the adopted 
neutron lifetime is t„ = 881.9 s [2J|. Both proton and 
neutron captures and other nuclear reactions of X-nuclei 
are taken into account. We have modified most of the 
thermonuclear reaction rates on the X-nuclei from the 
original rates (without X-nuclei). Binding energies be- 
tween nuclei and X particles are of order ~ 10 MeV. 
They are larger than those between a nucleus and mas- 
sive particles which only interact electromagnetically of 
~ 0.1 — 1 MeV [25l|. Hence, X particles lead to signifi- 
cant changes in reaction Q- values for reactions involving 
X-nuclei. Three possible effects of the binding of the 
X particles are: 1) changes in the Coulomb barriers re- 
sulting from the charge (if any) of the X particle in the 
nucleus; 2) modified reduced masses; and most impor- 
tantly, 3) the modified Q-values. 



A. Properties of the X Particle 

The X° particle is assumed to be hadronic and to have 
zero electric charge and zero spin. Its mass is assumed 
to be much larger than the nucleon mass. We note that 
X+ particles may also be present during BBN. Unlike 
the leptonic X+ case, they could have a strong interac- 
tion with nuclei. However, there exists Coulomb repul- 
sion leading to a certain degree of suppression of their 
reaction rates. Nevertheless, X+ particles should even- 
tually be included though we neglect them in the present 
investigation. Also, in this study, the nonthermal nucle- 
osynthesis triggered by the later electromagnetic and/or 
hadronic decay is not studied. These effects will be ad- 
dressed in a future publication. For now, however, as a 
first step we focus only on the effects of particles on 
BBN. 



B. Nuclear Binding Energies 

The nucleosynthesis of X-nuclei is strongly dependent 
upon their binding energies. In our calculations, bind- 
ing energies and eigenstate wave functions of X-nuclei 
are computed taking into account the nuclear interac- 
tion and Coulomb interaction between the nucleus and 
the X particle. We assume that the potential is spheri- 
cally symmetric. We then solve the two-body Shrodinger 
equation by a variational calculation (using the Gaussian 
expansion method \2d\) to obtain binding energies. 

The two-body Shrodinger equation for a spherically- 
symmetric system is 

(^-^S/^ + V{r)-E^ij{r)=0 , (3) 

where h is Planck's constant, fi is the reduced mass, V{r) 
is the central potential at r, E is the energy, and V'(^) 
is the wave function at r. Under the assumption that 
the X particle is much heavier than the light nuclides, 
fi is approximately given by the mass of the nuclide now 
considered. The central potential V{r) is composed of 
a nuclear interaction V^n('') and a Coulomb interaction 
Vb(r), i.e., 

V{r) = V^{r)+Vc{r) . (4) 

For the Coulomb potential, we assume that the charge 
distributions of the nuclei are Gaussian. We then use the 
charge radii determined from experiments of the corre- 
sponding nuclei (or neighboring nuclei when experimen- 
tal data are not available). Here, we write 

Vcir)^^^eri(^), (5) 

where Zac and Zxe are the charges of nuclide A and X, 
respectively. The parameter rp is related to the mean 
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square charge radius (r^) as rp = ■\/2/3(r^) ' , and 
erf(a;) = (2/0F) exp{-t^)dt. 

The X particle nucleon potential adopted here is as- 
sumed to be a square well of radius 2.5 fm and of depth 
of —25.5 MeV. This potential reproduces the binding en- 
ergy of the deuteron i.e., 2.224 MeV. A Woods-Saxon 
potential is adopted for the nuclear potential between 
other nuclei and the X-particles, i.e., 
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where the parameters are taken to be Vb = 50 MeV, 
a = 0.6 fm and R — (r^)^^^- The mean square matter 
radii for nuclei, (j'm)^ j taken from experiments of 
corresponding or neighboring nuclei. As a special case, 
the binding energy of two protons and an X particle sys- 
tem, i.e., ppxy in ^ Is orbit, is calculated with the same 
nuclear potential for the px and p system as that adopted 
for the estimation of the binding energy of and an X . 

The adopted radii and obtained binding energies are 
listed in Table [H Binding energies in the case of neu- 
tral X^, negatively-charged X~ and positively-charged 
X^ particles are shown in columns 6, 7 and 8, respec- 
tively. The adopted root mean square (RMS) nuclear 
matter radii and their references are listed in columns 
2 and 3. RMS charge radii and their references are 
shown in columns 4 and 5. Since the X particles are 
bound strongly to nuclei, their binding energies are typ- 
ically large (~ 10 MeV), and are even larger for heav- 
ier nuclei. Hence, they are bound to nuclei from early 
in the BBN epoch. We note that their binding ener- 
gies would be smaller 0.1 — 1 MeV) if they could 
only bind electromagnetically to nuclei. In that case 
they would not be bound to nuclei until low tempera- 
ture (Tg = T/(109 K) < 0.3). The obtained binding en- 
ergies are used for the estimation of Q- values of various 
reactions as described below. 



C. Reaction Rates 

1. Radiative X Capture Reactions 

We assume that the rates of radiative neutral X° cap- 
ture reactions by nuclei are roughly given by those of 
radiative neutron capture reactions by the nuclides or 
neighboring nuclides (if there are no corresponding data) . 
This assumption is introduced because we suppose that 
the X particles interact as strongly as normal nucleons. 
We correct the reduced mass and net charge for reactions 
involving X particles using the equations written below 
[Eqs. d?]), dH)]. The adopted reaction rates Na{(Jv), per 
second per mole cm~^, are shown in TableHH where Na is 
Avogadro's number. Reaction Q-values are derived tak- 
ing account of the binding energies of the X-nuclei listed 
in Table m 

There are two noteworthy cases, n{X ,"f)nx and 
p{X,j)px- In the n plus neutral A"" system, the electric 



multipole transitions do not occur because of charge neu- 
trality. The magnetic dipole transition also disappears 
by the orthogonality condition between the scattering- 
and bound-s-wave states. Although only the magnetic 
quadrupole or higher multipole transitions are allowed, 
they arc hindered by more than a factor of ~ 10^ com- 
pared with allowed electric dipole transition for a photon 
energy of a few MeV. In the n plus charged X^ system, 
the electric multipole transitions are allowed, but their 
transition probabilities disappear in the limit of a very 
massive X particle (m 3> 1 GeV). This is because A- 
multipole moment is proportional to m~^. The electric 
dipole transition rate, then, is very small for the reaction 
n{X ,"f)nx- Hence, we set the n(X ,"f)nx rate to zero. 

The nuclear potential for protons adopted in this study 
(Sec. IIIB|) leads to only one bound L — state with a 
binding energy of —9.2 MeV. Nuclei heavier than the 
nucleon can bind to the X particles in L = 1 excited 
states. In the system of p plus X, states exist with spin 
and parities of J'^ = 1/2+ (p) and 0+ {X), thus leading 
to a bound state with 1/2+. There is then no possibility 
for an electric dipole transition, i.e., spin change AL = 1 
and a parity change from a s-wave relative orbital angular 
momentum between a p and an X. The electric dipole 
transition to the bound state from a p-wave between the 
p and X is thus the dominant channel for the radiative 
capture reaction of a proton by an X particle. The rate 
of the p{X,^)px reaction has been estimated using the 
code RADCAP published by Bertulani [37*1 adopting the 
potential between a proton and an X particle as given in 
Secimi 



2. Nonresonant Neutron Capture Reactions of X -Nuclei 

We include reactions between neutrons and AT-nuclei in 
the reaction network. We adopt known reaction rates for 
normal nuclei whenever possible. When the correspond- 
ing reaction rates are not available, rates of reactions 
for neighboring nuclei are adopted. For neutron capture 
reactions at low energies, the s-wave interactions would 
dominate and the cross sections are proportional to the 
square of the de Broglie wavelength A = h/{iiv). We 
correct for the reduced mass and obtain neutron-induced 
reaction rates {av)Ax+n given by 



A 



(7) 



where {av)A+n is the neutron-capture reaction rate for 
the normal nucleus. A and Ax are the reduced masses 
for normal nuclei plus a neutron, and an A"-nucleus plus 
a neutron, respectively, in atomic mass units. 

In the case of radiative neutron capture, i.e., 
Ax{n,j)Bx , the electric dipole moment is very small 
similar to the n{X ^"f)nx reaction. Hence, the higher elec- 
tric quadrupole or magnetic dipole transitions contribute 
to the cross sections which are hindered by a factor of 
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TABLE I: Binding Energies of X Particles to Nuclei 



nuclide 


rS^^ (fm)" 


Ref. 


RMS /f \ /) 

rf ^ (fm) 


Ref. 


A case 


Ebm (MeV) 
X~ case 


X case 




- 


- 


0.875 ± 0.007 


[23 


9.242 


10.103 


8.391 




1.971 ± 0.005 


[28] 


2.116 ± 0.006 


[29] 


24.570 


25.344 


23.798 


TT 

■^Hx 


1.657 ± 0.097^= 




1.755 ± 0.086 


[30J 


24.013 


24.937 


23.091 




- 


- 


- 


- 


24.479 


26.312 


22.665 


^Hex 


1.775 ± 0.034'= 


M 


1.959 ± 0.030 


[30] 


25.819 


27.526 


24.117 


4xx 

Hex 


1.59 ± 0.04 


[31J 


1.80 ± 0.04 


[31] 


25.621 


27.491 


23.756 


^IXX 

'Hex 


2.52 ± O.OS*" 


[31J 


2.38 ± 0.03*= 


[31] 


38.221 


39.697 


36.748 


fixx 

'^Hex 


2.52 ± 0.03 


[31J 


2.38 ± 0.03 


[31] 


39.235 


40.724 


37.748 


^ X . 


2.35 ± 0.03^ 


[31J 


2.48 ± 0.03^ 


[31] 


36.763 


38.923 


34.607 


fix . 


2.35 ± 0.03 


[31J 


2.48 ± 0.03 


[31] 


37.853 


40.031 


35.679 


'^Lix 


2.35 ± 0.03 


[31J 


2.43 ± 0.02 


[31] 


38.695 


40.924 


36.470 


S X . 


2.38 ± 0.02 


[31J 


2.42 ± 0.02 


[31] 


39.610 


41.856 


37.368 


fix-. 

®Bex 


2.33 ± 0.02» 


[31J 


2.52 ± 0.02* 


[31] 


37.682 


40.551 


34.819 


'^Bex 


2.33 ± 0.02 


M 


2.52 ± 0.02 


[31] 


38.528 


41.415 


35.647 


*Bex 


2.33 ± 0.02» 


[31J 


2.52 ± 0.02* 


[31] 


39.203 


42.104 


36.307 


^Bex 


2.38 ± 0.01 


[31] 


2.50 ± 0.01 


[31] 


40.153 


43.080 


37.231 


^"Bex 


2.28 ± 0.02 


[31J 


2.40 ± 0.02 


[31] 


39.858 


42.912 


36.810 


^Bx 


2.45 ± O.IO** 


[32] 


2.68 ± 0.12'' 


[32] 


39.499 


42.910 


36.095 


*Bx 


2.45 ± 0.10 


[32] 


2.68 ± 0.12 


[32] 


40.138 


43.565 


36.717 


®Bx 


2.45 ± O.IO*" 


[32] 


2.68 ± 0.12'' 


[32] 


40.662 


44.102 


37.228 


^°Bx 


2.45 ± O.IO*" 


[32] 


2.68 ± 0.12'' 


[32] 


41.108 


44.560 


37.663 


Tlx-. 

"Bx 


2.45 ± O.IO*" 


[32] 


2.68 ± 0.12'' 


[32] 


41.490 


44.951 


38.034 


1 9 x-i 

^^Bx 


2.35 ± 0.02 


[31J 


2.51 ± 0.02 


[31] 


41.148 


44.835 


37.469 


10 

Ox 


2.32 ± 0.02' 


[31J 


2.51 ± 0.02' 


[31J 


40.170 


44.574 


35.777 


"Cx 


2.32 ± 0.02' 


[31J 


2.51 ± 0.02' 


[31] 


40.577 


44.995 


36.170 


19.^ 

^^Cx 


2.32 ± 0.02 


[31J 


2.51 ± 0.02 


[31] 


40.929 


45.359 


36.510 


^^Cx 


2.28 ± 0.04 


[33J 


2.463 ± 0.004 


[3^ 


40.955 


45.476 


36.445 


"Cx 


2.30 ± 0.07 


[33] 


2.496 ± 0.002 


[34] 


41.382 


45.857 


36.917 




2 47 + 07 


f33l 


2 62 + 07' 


f33l 


41.952 


46.906 


37.012 


"Nx 


2.31 ± 0.04 


[33j 


2.47 ± 0.04' 


[33] 


41.173 


46.430 


35.931 


"Nx 


2.47 ± 0.03 


[33J 


2.56 ± 0.01 


[3^ 


42.494 


47.572 


37.431 


I'Nx 


2.42 ± 0.10 


[33j 


2.61 ± 0.01 


[3^ 


42.420 


47.427 


37.425 


<Jx 


2.40 ± 0.03 


[33j 


2.56 ± 0.03' 


[33] 


42.058 


47.886 


36.247 


Ux 


2.44 ± 0.04 


[33J 


2.59 ± 0.04' 


[33] 


42.543 


48.302 


36.802 


Ux 


2.46 ± 0.12 


[35J 


2.71 ± 0.02 


[35] 


42.871 


48.412 


37.343 



"Root mean square (RMS) nuclear matter radius. 
''RMS charge radius. 

'^Derived by (r^^f = (rRMS)2 _ (aRMS)2 ^j^h qRMS = o.875± 
0.007 fm using a RMS proton matter radius determined in experi- 
ment as a RMS charge radius. 

''Estimated from the binding energies of ^Hx plus Q-values of the 
reaction 1Hx(p,7)pPx (See text in Sec. Ill Cj3) l. 

"^Taken from ^He radius. 

■'Taken from ^Li radius. 

^Taken from ^Be radius. 

''Taken from *B radius. 

'Taken from l-^C radius. 

'Derived by (rRMS)2 ^ (rRMS)2 + (aRMS)2 ^j^j, a^MS = 0.875 ± 
0.007 fm using a RMS matter radius determined in experiment. 
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TABLE 
A{X,^)A 



Rates of X'^ Radiative Capture Reactions 



Product 




Reaction Rate (cm s ^ mole 


Ref. 


Rev. Coef." 


Tlx 







0.987 


1 TT 

Hx 


4 X lO'' 


ft 


0.987 


2tt 

Hx 


7.4(1 + 18.9Tg) 




2.79 


^Hx 


4.2 X 10^ 


(St . 


5.13 


"^Hex 


4.1 X 10"^(1 + QOSTg) 


TT 


5.13 


Hex 


2.3 X 10^ 


(St • 


7.89 


(St • 

^Lix 


1.0 X 10^ 


flT • 


14.50 


Lix 


7.7 X 10^ 


7t . 


18.27 


*Lix 


5.9 X 10^ 




22.33 


^Bex 


7.6 X 10^ 


(St • 


18.27 


®Bex 


1.0 X 10^ 


^Be'' 


26.64 


*Bx 


5.9 X 10^ 


(St . 


22.33 


"Bx 


5.5 X 10^ 




31.20 


"Bx 


5.1 


1 1 T-» 

"B 


36.00 


Bx 


7.1 X 10"^ 


"c 


41.02 


"Cx 


4.5 X 10^ 


iUg 


36.00 


Cx 


2.7 


19^ 


41.02 


"Cx 


6.1 X 10"^ 


"C 


46.25 


"Cx 


5.2 X 10"^ 


"c 


51.69 


''Nx 


3.8 X 10^ 


lOg 


41.02 


"Nx 


2.3 


12^ 


46.25 


"Nx 


4.4 X 10^ 


"N 


51.69 


i^Nx 


1.3 X 10"^ 




57.32 


"Ox 


2.0 




51.69 


Ux 


3.8 X 10^ 


"N 


57.32 




8.7 X 10-2 


16qc 


63.15 



"For nucleus i with mass number Ai, the reverse coefficient is 
defined as in [36ll . They are given by 0.9867A^''2 for the process 
i(X,^)ix on the assumption that the X particle is much heavier 
than a nucleus. 

'Approximate values calculated with a code RADCAP *37*| at 
temperatures Tg ~ 2 — 6. 
<^ Taken from Ref. [H. 



^ 10'^ for emitted photon energies of order ^ lOMeV. 
We adopt corresponding reaction rates for normal nu- 
clei multiplied by 10^^ to account for this hindrance 
of radiative capture cross sections. Especially, the rate 
of the ^Hx(f^,7)^Hx reaction becomes negligibly small 
compared with that of the ^Y{x{n,p)nx reaction (see Ta- 
ble |TTl|), which predominantly processes ^Hx- 



3. Nonresonant Reactions Between Charged Particles 



The leading term in the expression for thermonuclear 
reaction rates (TRR) {av) between charged particles can 



be roughly written (e.g., [33]) as 

7.1/2 



MNK^(^)^^^-^.(i.„)exp(-.) , 



(8) 



where Eq = 1.22{zfZ^AT^y/^ keV is the energy at the 
peak of the Gamow window, Afu is the atomic mass unit, 
S{Eo) is the "astrophysical S'-factor" at Eq, is the 
Boltzmann constant, Tg is the temperature in units of 
10^ K, and 



42.46 



zJZlA 



1/3 



(9) 



The astrophysical S factor contains the nuclear matrix 
element for the reaction. We assume that the S{Eo) 
values for reactions involving X-nuclei are the same as 
those for the reactions of the corresponding normal nu- 
clei [1^ |4l[ . When the corresponding S factors are 
not available, S factors of reactions for neighboring nuclei 
are adopted. Corrections for the TRR in the above equa- 
tion arise from the reduced mass A. z\ and Z2 are the 
atomic numbers for the projectile normal nucleus and the 
target X-nucleus, respectively. Note, that we have taken 
the spin of the X particles to be zero in this study. 

Rates for two reactions, px{p,7)ppx and nx(p,7)^Hx, 
are calculated with the code RADCAP [s^l ■ The adopted 
nuclear potential is that of the d plus X system (Sec. Ill B"|) 
for both reactions. 

The adopted reaction rates iVyi(o't;)NR, thus obtained, 
are shown in Table IIIII for radiative reactions. Table IIVI 
is for nonradiative reactions independent of the deuteron 
and Table |V] is for deuteron capture nonradiative reac- 
tions. 



4- Rates for Reactions with Negative Q-Values 

We find that there are several reactions whose Q- 
values become negative when nuclei are bound to X par- 
ticles. We neglect most of those reactions except for the 
3Hx(p,n)3Hex, ^llex{t,nf Lix and ^Bex{^'Re,p)^Ux 
reactions. However neglected reactions might be impor- 
tant and should eventually be included in our calcula- 
tion. The rates of these three reactions are given by 
the Hauser-Feshbach approximation as follows. For a 
compound-nucleus reaction 



1 



C 



(10) 



where C is the compound nucleus, the cross section is 
given by the product of the probability of formation of 
the compound nucleus from 1-1-2, and that of its decay 
into the 3-1-4 particle channel, i.e. 



(1 + (5i2)7rA 



1 



i2(2/i + l)(2/2 + l) 
X Yl l(3,4|i/„|C)(C|i7,|l,2)r 



(11) 
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TABLE III: Radiative Reaction Rates for X°-Nuclei 



Reaction 


Reaction Rate (cm^ s ^ mole ^) 


Reverse Coefficient'' 


Q (MeV) 


^nx{p,'y)^1ix 


4 X 10^ 


1.32 


17.553 


iHx(n,7)'Hx 


1.2 X 10^ 


1.32 


17.553 


^Hx(p,7)^PPx 


8 X 10'^r9"^/^exp(-4.25/rg''^)' 


3.95 


15.237 


2Hx(n,7)^Hx 


2.9 X 10"^(1 + IS.QTg) 


2.96 


5.700 


'Hx(p,7)'Hex 


2.3 X 10^r9"^''^exp(-4.25/rg'^^) 


2.96 


6.742 


='Hx(p,7)^Hex 


2.0 X lO^Tg'^^'^ exp{~4.25/T^^'^) 


3.95 


21.422 


^Hex(n,7)''Hex 


3.7 X 10"^(1 + 905.r9) 


3.95 


20.379 


^Hex(a,7)^Bex 


3.6 X 10''rg"^/^exp(-16.99/Tg''^) 


3.95 


17.346 


^Hex(n,7)5Hex 


3.7"^ 


0.493 


8.656 


*Hex(p,7)''Lix 


5.6 X 10'rg"'''^exp(-6.74/r9/^)'' 


0.493 


9.176 


'*Hex(d,7)'^Lix 


2.6 X 10^r9"^''^exp(-8.50/rg'^^) 


2.79 


13.706 


''Hex(t,7)'Lix 


2.5 X 10^Tg"^''^exp(-9.73/rg'^^) 


2.56 


15.541 


*Hex('He,7)"Bex 


4.0 X 10^r9"^/^exp(-15.44/r9/^) 


2.56 


14.494 


*Hex(a,7)*Bex 


3.6 X 10^Tg"^'^%xp(-16.99/rg''^)' 


7.89 


13.490 


*Hex('Li,7)iOBx 


3.0 X 10^rg"^/^exp(-25.49/r9''^) 


6.22 


19.949 


5Hex(n,7)*'Hex 


3.7^= 


7.89 


2.880 




5.6 X 10'5Tg"^^^exp(-6.74/r9''^)'^ 


2.63 


4.224 


'5Hex(p,7)'Lix 


5.6 X 10'^rg"^^^exp(-6.74/r9/^)'^ 


0.493 


9.436 


''Uxin,jfUx 


3.7^^ 


2.63 


6.754 


'^Lix(p,7)''Bex 


6.4 X 10''rg"'/^exp(-8.84/r9''^)'' 


7.89 


1.513 


«Lix(n,7)^Lix 


3.7 


1.48 


8.092 


''Lix(p,7)^Bex 


6.4 X 10^rg"^/^exp(-8.84/rg''^) 


1.48 


6.281 


«'Lix(a,7)"'Bx 


3.4 X 10''Tg"^/%xp(-22.27/Tg/^) 


3.38 


7.716 


^Ux{n,"ffUx 


3.8 


1.58 


2.947 


^Lix(p,7)*Bex 


1.7 X 10^rg"^/^exp(-8.84/rg''^) 


7.89 


17.763 


"Lix(a,7)"Bx 


3.1 X 10'^rg"^''^exp(-22.27/rg''^) 


7.89 


11.460 


*Lix(p,7)^Bex 


1.7 X 10'^rg"^/^exp(-8.84/r9/^/ 


2.47 


17.431 


«Bex(n,7)^Bex 


3.7'= 


0.493 


11.522 


^Bex(n,7)»Bex 


3.7"= 


7.89 


19.574 


^Bex(p,7)*Bx 


3.0 X 10^r9'^/^exp(-10.71/r9''^) 


1.58 


1.747 


j3ex(,Q!,7j 


1 .O X lU Ig ^Xp^^ A0.OI/1q ) 




y.ooy 


*Bex(n,7)''Bex 


3.8" 


0.493 


2.615 


*Bex(p,7)^Bx 


3.0 X 10^rg"^''^exp(-10.71/rg''^)'' 


0.493 


1.274 


''Bex(n,7)i«Bex 


8.2 X 10"^' 


7.89 


6.518 


SBex(p,7)'''Bx 


1.3 X 10'^rg"^/^exp(-10.71/r9/^) 


1.13 


7.542 


i''Bex(p,7)"Bx 


1.3 X 10^rg~^''^exp(-10.71/r9''^)^ 


0.493 


12.859 


"For nuclides a = k, .. 


. with mass numbers Aa and numbers of 







magnetic substates , the reverse coefficients are defined as in |36|I : 
on the assumption that an X particle is much heavier than nuclides, 
they are given by 0.9867{gigj / gk)A'^^^ for the process ixij ,'y)kx ■ 

'The approximate vales calculated with a code RADCAP at 
temperatures Tg ~ 2 — 6. 

■^The rate of the reaction ®Li(n,7)^Li multiplied by 10"'' is used. 

"^The S-factor for the reaction ^Li(p,7)'^Be is used. 

"^The S-factor for the reaction ^He(a,7)^Be is used. 

■'^The S-factor for the reaction '^Li(p,7)®Be is used. 

"The rate of the reaction '^Li(n,7)*Li multiplied by 10"'' is used. 

''The S-factor for the reaction '^Be(p,7)*B is used. 

'The cross section from l38l multiplied by 10~^ is used. 

^The S-factor for the reaction ^Be(p,7)^''B is used. 
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TABLE III: Continued. 



Reaction 


Reaction Rate (cm"^ s ^ mole ^) 


Reverse Coefficient" 


Q (MeV) 


«Bx(n,7)«Bx 


3. 7* 


2.47 


19.101 


^Bx(n,7)'°Bx 


3.8" 


1.13 


8.883 




4.5 X 10^T~^^^ ex.p{-12A2/Tg^'^y 


7.89 


3.514 


"Bx(n,7)"Bx 


5.5 X 10^ 


3.45 


11.835 


"Bx(p,7)"Cx 


4.5 X 10^r3"^/^exp(-12.42/rg''^) 


3.45 


8.158 


"Bx(n,7)i2Bx 


6.1 X 10"^ 


2.63 


3.029 


"Bx(p,7)''Cx 


4.5 X lO'^Tg^^'-^ exp{~12A2/TQ^^) 


7.89 


15.396 


i°Cx(n,7)"Cx 


3.8™ 


0.493 


13.527 




5.5 X 10^" 


7.89 


19.074 


"Cx(p,7)''Nx 


4.1 X 10*r3"^''^exp(-14.03/rg''^) 


2.63 


1.977 


i2Cx(n,7)''Cx 


3.8 X 10"^ 


0.987 


4.972 


^2Cx(p,7)"Nx 


2.0 X 10^r3"^''^exp(-14.03/ry'^^) 


0.987 


2.187 


i2Cx(a,7)''Ox 


9.4 X 10^rg"^/^exp(-35.35/rg''^) 


7.89 


9.104 


"Cx(n,7)"Cx 


1.0 X IQ-i 


3.95 


8.603 


"Cx(p,7) Nx 


7.8 X lO^Tg ^•^exp(-14.03/rg'"^) 


1.32 


9.089 




fi R V 1 n^T"'^/'^ pYnf' ^An'^/T^^^^ 

u.u A ±u -t g exp^ — L^.yjOf ± g j 


u. yo * 


1 1 94-"^ 


"Nx(p,7)"Ox 


3.9 X 10^rg"^/^exp(-15.55/Tg^^) 


3.95 


5.512 


"Nx(n,7)i5Nx 


8.7 


2.96 


10.759 


"Nx(p,7)''Ox 


4.8 X 10^rg"^/^exp(-15.55/rg''^) 


2.96 


7.346 


15Nx(p,7)"Ox 


9.6 X 10**rg"^/^exp(-15.55/rg''^) 


3.95 


12.578 


"For nuclides a = i,j, k, .. 


. with mass numbers Aa and numbers of 







magnetic substates Qa , the reverse coefficients are defined as in |36|I : 
on the assumption that an X particle is much heavier than nuclides, 
they are given by 0.9867(gigj /gk)A^y^ for the process «x(ii7)fcx. 

*The rate of the reaction '^Li(n,7)®Li multiplied by 10"'^ is used. 

*The cross section from [38l | multiplied by 10~^ is used. 

'The 5-factor for the reaction ^"6(^,7)^^0 is used. 
"The rate of the reaction '^Li(n,7)*Li multiplied by 10^'^ is used. 

"The rate of the reaction ^''B(n,7)^^B multiplied by 10~^ is used. 



where li is the spin of the nucleus i, and A12 is the 
de Broghe wavelength of the entrance channel (e.g., Eq. 
(3.9.26) in Ref. [H) satisfying 
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Ai-jE. 



barn. 



(12) 



The factor (1 + Sij) doubles the cross section for in- 
distinguishable particles. The matrix elements in angle 
brackets have information on the nuclear factors (and 
the Coulomb barrier penetration probabilities if they are 
reactions between charged particles). The barrier pene- 
tration probability for s- waves in the low energy limit is 
given (e.g., Eq. (3.10.10) in Ref. [H) by 



exp I h 4 



Er 



1/2N 



1/2 



exp -31.28- 



ZiZiA 



1/2 



E, 



1/2 



1 /2 

+1.05 (Aij^xRij.fmZtZj) 



(13) 



where ZiC and Zje are the charges of nuclei i and jx, 
respectively, Aij^x is the reduced mass, Eij and i^i j,kcV 
denotes the center of mass energy, where units of keV 
are indicated where used. Vij is the relative velocity 
of the projectile target system of i and jx- Rij — 
lA^Ay^^ + A^^) fm is the contact nuclear radius, i.e. 
the separation between the centers of particle i and jx 
when the attractive nuclear force overcomes the Coulomb 
barrier. Here Ai and Aj are the masses of nuclides i and 
j in atomic mass units. Rij^im is the Rij value in units 
of fm. Ec = 1.44 ZiZj/Rij^tra MeV is the height of the 
Coulomb barrier. The penetration probabilities are con- 
tained in the matrix elements. Defining partial widths Ta 
and Fb of the compound nucleus for decays into entrance 
and exit channels, respectively, the cross section for the 
reaction Eq. ITUl) has a scaling relation of 



CTj, cx A,,rarb cx 



TgTb 

A- F- 



(14) 
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TABLE IV: Nonradiative Reaction Rates for X"-Nuclei 



Reaction 


Reaction Rate (cm s mole ) 


Reverse Coemcient 


Q (MeV) 


Itt / \1 

iix(n,p) nx 


2.0 X 10^' 


1.0 


0.0 


2 / \ 2tt 

ppx(n,p) Hx 


2.0 X 10^' 


0.333 


2.316 


^'Hx{a,pfUex 


2.8 X lO^^Tg"^'''' exp( — lO.Tl/Tg''^)'' 


8.0 


11.150 


^Hx(p,n)''Hex 


2.5 X lO^^Tg"^/^ exp(-4.25/rg^'^^) 


1.0 


1.043 


^}iexit,nfLix 


8.7 X 10^''rg"^/^exp(-9.73/Tg^''^) 


1.73 


7.449 


4tt /3tt \6t • 

*Hex('^He,p)''Lix 


1.1 X lO-'^-'^Tg"^^^ exp( — 15.44/rg''^) 


1.73 


8.213 


5 T • / \ n T T 

''Lix {n,pymex 


2.0 X 10^' 


1.0 


2.530 


^Lix(p,Q)''Hex 


1.0 X 10''Tg"^/^exp(-8.84/r9''^) 


1.00 


4.273 


''Ux{p,nyBex 


8.3 X lO^Tg"^^^ exp( — ll.lS/Tg''^)'^ 


5.0 


14.816 


8Lix(p,Q)^Hex 


1.0 X lO^Tg"^''^ exp(-11.13/rg^'^^)'' 


0.313 


13.032 


8t • / \ 11 T-i 

'*Lix{a,n)"^Bx 


7.5 X lO^^Tg"^''^ exp(-22.27/rg''^) 


5.0 


8.512 


^Bexin,pfUx 


2.0 X lO'" 


0.333 


5.241 


'^Bex(n,p)^Lix 


2.0 X 10^ 


1.00 


1.810 


"Bex(a,n)^^Cx 


4.1 X 10^^rg"^''^exp(-26.98/rg''^) 


16.00 


6.477 


*'Bx(n,p) Bex 


3.2 X 10*^ 


5.0 


17.827 


''BxCa.p) Cx 


9.4 X lO^^Tg"^''^ exp(-31.30/rg''^) 


5.00 


7.846 


^Bx(n,p)^Bex 


2.0 X 10^* 


1.0 


1.341 


^Bx(n,Q:)^Lix 


4.2 X 10*^ 


0.333 


13.539 


^"Bx(a,n)"Nx 


1.1 X lO^^Tg"^''^ exp(-31.30/rg''^) 


14.0 


1.123 


lO-i-i / \13,^ 

^"Bx (q;,p)"Cx 


8.6 X 10"rg"^''^exp(-31.30/rg''^) 


14.0 


3.908 


^'Bx(p,a)'*Bex 


2.0 X lO^^Tg"^''^ exp(-12.42/rg^''^)'' 


1.0 


6.303 


^^Bx(a,n) Nx 


6.3 X lO^^Tg"^^^ exp(-31.30/rg^^) 


5.33 


1.162 


"Bx(a,p) Cx 


4.8 X 10"rg"^''^exp(-31.30/rg''^) 


16.0 


0.676 


^^Bx(p,n)'^Cx 


3.9 X lO^^Tg"^^^ exp( — 12.42/rg''^) 


3.00 


12.367 


^^Bx(p,a)-^Bex 


2.0 X 10"rg"^/^exp(-12.42/rg''^) 


0.188 


5.889 


^^Bx(a,n)^''Nx 


2.8 X 10^^rg"^''^exp(-31.30/r9^^) 


6.00 


8.892 


^"Cx(n,p)^"Bx 


2.0 X 10^' 


0.143 


5.369 


^"Cx(n,a)^Bex 


4.2 X 10**^ 


0.0625 


3.933 


"Cx(n,p)"Bx 


1.4 X 10* 


1.0 


3.677 


11 / \8i-. 

"Cx(n,a)'*Bex 


1.3 X 10*' 


1.0 


9.981 


"Cx(a,p) Nx 


6.4 X 10^''^r9"^^^exp(-35.35/r9''^) 


5.33 


4.840 


"Cx(a,n)^''Ox 


6.2 X 10^^r9"^''^exp(-35.35/r9^^) 


8.00 


4.131 


12at / \12 

^^Nx(n,p)'^Cx 


1.4 X 10*^ 


3.0 


17.097 


^^Nx(n,a)"Bx 


4.2 X 10** 


0.188 


9.278 


^^NxfQ,p)^^Ox 


5.1 X 10^''r9"^''^exp(-39.18/r9^^) 


6.00 


10.209 


i3Nx(n,p)"Cx 


1.6 X 10* 


1.00 


2.049 


"Nx(a,p)^<^Ox 


3.0 X 10^'^r9"^''^exp(-39.18/r9''^) 


8.00 


6.916 


i5Nx(p,a)''Cx 


1.1 X 10^^r9"^/^exp(-15.55/Tj/^) 


0.500 


3.474 


i50x(n,p)''Nx 


3.1 X 10* 


1.00 


3.413 


^50x(n,a)i2Cx 


3.1 X 10^ 


0.500 


6.888 



"For nuclides a = i,j,k,... with mass numbers Aa and numbers 
of magnetic substates ga, the reverse coefficients are defined as 
in [ssll : on the assumption that an X particle is much heavier than 
nuclides, they are given by igigj/{gk9l)){Aj l-^h)^^'^ for the process 
«x(i,fc)Zx- 

'The rate of the reaction ^Be(n,p)^Li is used. 
'^The S-factor for the reaction ®B(ci,p)^^C is used. 
■^The S-factor for the reaction *Li(p,na)'*He is used. 
'^The 5-factor for the reaction ^Li(p,o)'*He is used. 
■''The rate of the reaction *B(n,pa)'*He is used. 
SThe rate of the reaction ^''Be(n,o)^Li is used. 
''The S-factor for the reaction ^^B(p,a)®Be is used. 
'The rate of the reaction ^^C(n,2a)*He is used. 
^The rate of the reaction ^^C(n,p)'^^B is used. 
*The rate of the reaction ^''B(n,a)'^Li is used. 
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TABLE V: (d,n) and Reaction Rates for X°-Nuclei 



Reaction 


Reaction Rate (cm'^ s ^ mole ^) 


Reverse Coefficient" 


Q (MeV) 




3.3 X 10*rg"^''^exp(-5.35/rg''^) 


6.364 


2.182 


2Hx(d,p)'Hex 


3.1 X 10**rg"^''^exp(-5.35/rg''^) 


6.364 


5.811 


3Hx(rf,n)4Hex 


9.0 X 10^"r9~^/^exp(-5.35/rg/^) 


8.485 


19.197 


^Rex{d,pfRex 


4.2 X 10^"r9"^''^exp(-8.50/rg/^) 


8.49 


18.154 


^Rex {d,nf Lix 


1.1 X 10"rg"^''^exp(-8.50/rg''^)'' 


1.061 


6.952 


''Rex{d,pfRex 


4.2 X 10"r9~^''%xp(-8.50/rg''^)'= 


1.061 


9.482 


'^Rex{d,nfLix 


1.1 X 10"r9"^''^exp(-8.50/rg''^)'' 


5.657 


2.000 


''Rexid,pfRex 


4.2 X 10"r9"^''^exp(--8.50/rg^^)'= 


16.97 


0.655 


'^Rex{d,nyUx 


2.3 X 10^^r9"^''^exp(-8.50/r9^^)'' 


1.061 


7.211 


''Ux{d,pfUx 


4.8 X 10^''r9"^''^exp(-11.13/r9^^)" 


5.657 


4.530 


^Ux{d,nyBex 


2.7 X 10"Tg~^/^exp(-11.13/rg''^)'' 


3.182 


4.056 


''Ux{d,pyUx 


8.9 X 10"r9"^''^exp(-11.13/r9''^)'= 


3.182 


5.867 


''Ux{d,nfBex 


2.7 X 10"Tg"^/^exp(-11.13/r9''^)'' 


16.97 


15.539 


''Ux{d,pfUx 


8.9 X 10"r9"^''^exp(-11.13/r9''^)'= 


3.394 


0.723 


^Ux{d,nfBex 


2.7 X 10"Tg"^/^exp(-11.13/r9''^)'' 


5.303 


15.206 


^Bex(rf,p)^Bex 


9.8 X 10"r9"^''%xp(-13.49/r9^^)'= 


1.061 


9.298 


''Bex{d,pfBex 


9.8 X 10"r9"^''%xp(-13.49/r9''^)'= 


16.97 


17.350 


''Bex{d,pfBex 


9.8 X 10"r9"^''%xp(-13.49/r9''^)'= 


1.060 


0.391 


^Bex(d,n)i»Bx 


4.2 X 10^^r9"^''^exp(-13.49/r9''^) 


2.424 


5.317 


^Bx{d,pfBx 


1.1 X io"r9"^''%xp(-i5.65/rg/^) 


5.303 


16.877 


»Bx(d,n)i»Cx 


4.5 X 10"Tg"^/^exp(-15.65/r9''^/ 


16.97 


1.290 


SBx(d,p)i''Bx 


1.3 X 10i2r9~^''^exp(-15.65/r9''^)3 


2.424 


6.658 


i°Bx(d,n)"Cx 


2.2 X 10^^Tg~^/^exp(-15.65/rg''^)'' 


7.425 


5.933 


i°Bx(d,p)"Bx 


1.3 X 10^^r9"^''%xp(-15.65/rg''^) 


7.425 


9.611 


"Bx(d,n)i2Cx 


2.2 X 10^^Tg"^/^exp(-15.65/rg''^)'' 


16.97 


13.172 


"Bx(d,p)''Bx 


1.3 X 10^^rg"^''^exp(-15.65/r9''^)» 


5.657 


0.805 


i"Bx(d,n)"Cx 


2.2 X 10^^Tg"^''^exp(-15.65/rg''^)'' 


6.364 


15.115 


i°Cx(d,p)"Cx 


1.4 X 10^^r9"^''^exp(-17.67/r9''^)» 


1.061 


11.302 


"Cx(d,p)''Cx 


1.4 X 10^^r9"^^%xp(-17.67/r9''^)» 


16.97 


16.849 




1.4 X 10^^r9"^''%xp(-17.67/r9''^)'' 


2.121 


2.748 


"Cx(d,n)"Nx 


2.3 X 10^^Tg"^/^exp(-17.67/rg''^)'' 


2.828 


6.865 




1.4 X 10^^r9"^''%xp(-17.67/r9''^)s 


8.485 


6.379 




2.3 X 10^^Tg"^/^exp(-17.67/r9''^)'' 


2.121 


9.021 


^^Nx(d,p) ■^Nx 


1.5 X lO^^Tg"^''^ exp(-19.59/rg''^)« 


6.364 


17.060 


i3Nx(d,n)"Ox 


2.4 X 10^^Tg"^/^exp(-19.59/r9''^)'' 


8.485 


3.288 


i3Nx(d,p)i*Nx 


1.5 X 10^^r9"^''^exp(-19.59/r9''^)'' 


2.828 


9.650 


"Nx(d,n)i5Q^ 


2.4 X 10^^Tg"^/^exp(-19.59/rg''^)'' 


6.364 


5.122 


i*Nx(d,p)i5Nx 


1.5 X 10^^r9"^''%xp(-19.59/r9^^)» 


6.364 


8.535 


i5Nx(d,n)«Ox 


2.4 X 10^^Tg~^/^exp(-19.59/rg''^)'' 


8.485 


10.354 



"For nuclides a = i,j,k,... with mass numbers Aa and numbers 
of magnetic substates ga, the reverse coefficients are defined as 
in [3^ : on the assumption that an X particle is much heavier than 
nuclides, they are given by {gigj / {9k9l ))(Aj /^k)^^"^ for the process 
«x(j,fc)'x. 

''The S-factor for the reaction ^H(d,n)^He is used. 
'^The S-factor for the reaction ^He(d,p)''He is used. 
"^The S-factor for the reaction '^Li(d,na)'*He is used. 
"^The S-factor for the reaction '^Be(d,pa)*He is used. 
■'^The S-factor for the reaction ^Be(d,n,)'^''B is used. 
*The S-factor for the reaction ^''Be(ii,p)^'^B is used. 
''The S-factor for the reaction ^^B{d,ny^C is used. 
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The partial width for the particle decay channel can be 
written [421 as 



(15) 



where 0^ is the dimensionless reduced width, i.e., a mea- 
sure of the degree to which the compound nuclear state 
can be described by the relative motion of i and j in a 
potential. For reactions of X-nuclei, we assume that the 
nuclear radius R and the reduced width 6^ are the same 
as those for the corresponding normal reactions. The 
cross section, therefore, scales according to: 



0-12 oc 



(wi2^'i2) (W34-P34) 



^12£'12 



(16) 



We use this scaling relation and adopt coefficients in Eq. 
()16|) from the standard nuclear reactions assuming that 
the coefficients contain the information of the purely nu- 
clear part and that other parts including Coulomb pene- 
tration factors related to corrected reaction Q-values can 
be extracted as in Eq. pB]) . 

The dimensionless reduced width 9^ is also related to 
the spectroscopic factor for a direct reaction [s^]. The 
distinction between a compound-nucleus and a direct re- 
actions can be obscure as low-energy direct reactions can 
result from many overlapping resonances. 

For the reaction '^Hx(p,n)'^Hex, we adopt the non- 
resonant rate of the normal reaction ■^He(n,»)^H, i.e., 
Na{<Jv)sbbn = 7.21 X 10* cm3 s'^ mole"! [Ilj. The 
penetration factor for the exit channel for ■^He(n,p)'^H 
is assumed to be Pp = 1 because of the high Q-value 
(Q > Ec)- Eq. leads to the following S'-factor for 
the '^Ilx{p,n)^llex reaction 



exp(— 27r?7) 



= 3.1 MeV barn, 



(17) 



where 77 = z\Z2e^ jijiv). 

For the reaction ^Hej(:(i,'T-)^Lix, we adopt the nonreso- 
nant rate of the ^Li(n,a)^H reaction, i.e., iVA(cr?;)sBBN — 
1.68 X 10* cm^ s"^ mole"^ [2l|. The penetration fac- 
tor of the exit channel for ^Li(n,a)'^H is also assumed 
to be Pq = 1 because of the high Q-value (Q > Ec)- 
Similarly, the S'-factor of the '*He;is:(t,n)^Lix reaction is 
derived from Eq. ([16]) to be 



cx+t 



11 MeV barn. 



(18) 



For the ^Hej(:('^He,p)^Lix reaction we adopt the non- 
resonant part of the S'-factor from the ^Li(p,Q;)'^He cross 
section, i.e., S'sbbn = 3.14 MeV barn [2l[. The penetra- 
tion factors of the exit channel for the ^Li(n,Q;)'^H and 
^Hex('^He,p)^Lijf reactions are assumed to be P^ = 1 
and Pp = 1 because of the high Q-values (Q > Ec)- The 
S'-factor for the ^Hejc ('^He,p)^Lix reaction is derived us- 
ing Eq. to be 



'S'-'Hex+^Hc = 63 MeV barn. 



(19) 



5. Transfer Reactions px (n,p)nx and px (ct,p)'^llex 

Since neutron radiative X capture reactions would be 
relatively weak, the most important reaction for neu- 
trons to become bound to X particles is px{n,p)nx- In 
this reaction, an X particle transfers from a proton to 
a neutron. For this reaction, we use the rate for the 
^Be(n,p)^Li reaction, which is similar to the px{n,p)nx 
in the sense that both ^Li and X are massive and strongly 
interacting spectator particles so that their reactions 
have similar dynamics. 

If Px were to survive beyond the epoch of *IIe produc- 
tion in the standard BBN, i.e., to temperatures Tg < 0.1 
(although this is found not to be the case in the present 
network calculation,) then ^He could become bound to 
an X particle via the reaction px{ct,p)'^ilsx- This kind 
of exchange reaction plays an important role in the cat- 
alyzed BBN scenario with only electromagnetically inter- 
acting X^ particles For the rates of this reaction, 
we use the ^B{a,p)^^C reaction since ''Be and X are mas- 
sive and strongly interacting spectators in the reactions. 
Furthermore the cross section is corrected for Coulomb 
penetration factors using Eq. The following rela- 
tion is then derived. 



~px+a 



5.1 X lO^^Ss 



i+a- 



(20) 



Based upon the S-factor for the *B(q;,p)^^C reaction, 
(S'sbbn = 8.88 x 10^ MeV barn,) we obtain Sp-^^a — 
45 MeV barn. 



6. P -Decay of X -Nuclei 

When Q- values are larger than the electron mass in /3- 
decay reactions, decay rates F scale as the fifth power of 
the Q- values, i.e., T oc Q^- To estimate the /3-decay rates 
for X-nuclei, we use /?-decay rates F of corresponding 
normal nuclei corrected for the phase-space factors, i.e., 



X 



Q 



In 2 

Tl/2 



Q 



(21) 



where is the half life of the normal nuclide, Qx and 
Q are the Q-values for the /?-decay of the X- and normal 
nuclide, respectively. We show the adopted /?-decay rates 
in Table IVll The ^IIe(/3~)^Li reaction rate is used to 
estimate the ^Bex{P'^)^^^x reaction. We neglect the re- 
actions ^UxiP+)^llex, ^Bex{f3+)'Ux, ^Bx{f3+)^Bex, 
and in the X+ case, px{f3~^)nx since the Q-values for 
these reactions are relatively small and their lifetimes 
would not be short compared to the BBN timescale. 

When the Q-value is < 1 MeV, atoms can decay pre- 
dominantly by electron capture [39|. However, the elec- 
tron capture can be neglected since we are considering 
only the high energy epoch of the early universe when 
the nuclei and X-nuclei are fully ionized. 



11 



TABLE VI: /3-Decay Rates of X-Nuclei 



Reaction 


X case 


Qx (MeV) 
X~ case 


case 


Q (MeV) 


2^1/2 


Ref. 


X case 


Decay Rate (s 
X~ case 


X"*" case 




^n(/3")^H 


0.782 


1.643 


-0.069 


0.782 


10.19 m 


[24J 


1.133x10" 


-3 


4.634x10"^ 


- 






1.825 


2.608 


1.044 


0.019 


(12.3±0.1) y 


M 


1.626x10 


1 


9.700x10^ 


9.954x10" 






0.726 


0.042 


1.409 


-0.732 


- 


[45] 


- 




- 


- 




®He(/3-)'^Li 


2.127 


2.815 


1.439 


3.508 


(806.7±1.5) ms 


m 


7.033x10" 


-2 


2.860x10"^ 


9.983x10" 


3 


®Be(/3+)'^Li 


3.437 


2.746 


4.126 


3.266 


- 


[45] 


7.757x10" 


-1 


2.525x10"^ 


1.934 




^Be(/3+)^Li 


0.006 


-0.651 


0.662 


-0.160 


- 


[45] 


- 




- 


- 




*Li(^-)*Be 


15.598 


16.253 


14.945 


16.005 


(839.9±0.9) ms 


[47] 


7.255x10" 


-1 


8.912x10"^ 


5.857x10" 




*B(/3+)*Be 


16.023 


15.497 


16.547 


17.980 


(770±3) ms 


[47] 


5.059x10" 


-1 


4.281x10"^ 


5.943x10" 




®B(/3+)''Be 


-0.463 


-0.976 


0.049 


0.046 


- 


[45] 


- 




- 


- 




1 n -M-v /■ -1 \ 1 n -M-i 

^°Be(/3 y°B 


1.806 


2.204 


1.409 


0.556 


(1.51±0.04)xl0^ y 


[47] 


5.263x10" 


12 


1.423x10"" 


1.520x10" 


12 


"C(/3+) "B 


3.564 


2.612 


4.512 


2.626 


(19.290±0.012) s 


[47] 


1.655x10" 




3.497x10"^ 


5.381x10" 




"C(/?+)^^B 


1.873 


0.917 


2.825 


0.960 


(1223. 1±1. 2) s 


[48] 


1.602x10" 


-2 


4.501x10"" 


1.250x10" 




"Bffl-V^C 


13.149 


13.892 


12.410 


13.370 


^20 20-1-0 02'! ms 


f48l 


3 157x10 


1 


4.155x10^ 


2 363x10 


1 


i2n(/3+)12c 


15.293 


14.768 


15.814 


16.316 


(11.000±0.016) ms 


[48] 


4.557x10 


1 


3.828x10^ 


5.389x10 


1 


"N(/3+)"c 


0.981 


0.245 


1.712 


1.199 


(9.965±0.0004) m 


[34] 


4.250x10" 


^4 


4.127x10"^ 


6.896x10" 


3 


"C(/3-)i*N 


1.268 


1.871 


0.670 


0.156 


(5730±40) y 


[34] 


1.340x10" 


-7 


9.377x10"^ 


5.519x10" 


9 


"0(/3+)"N 


4.558 


3.807 


5.305 


4.121 


(70.606±0.018) s 


[34] 


1.624x10" 


-2 


6.606x10"-^ 


3.470x10" 


2 




1.609 


0.857 


2.355 


1.732 


(122.24±0.16) s 


m 


3.922x10" 


-3 


1.683x10"** 


2.638x10" 


2 



D. Reaction Network 

The reaction network for bound X-nuclei is shown in 
Fig. [TJ Sohd arrows show nuclear reactions in the di- 
rection of positive Q-value while dashed arrows indicate 
f3^ decays. The network code includes reactions up to 
oxygen isotopes. However, we did not find significant 
nuclear flow beyond the nitrogen isotopes. Hence, this 
network code is more than large enough to calculate the 
evolution of the nuclear abundances. The adopted nu- 
clear reaction rates are summarized in Tables IHIl IIVI IVl 
and IVII In our code, radiative X-capture reactions are 
also included although they are not explicitly shown on 
Fig. [11 



III. RESULTS 

A. BBN Calculation Result 

Figures [2^ and show results of a BBN calculation 
for the case Yx = Nx/nt, — 10~^, where Nx and ni, are 
the number densities of the particles and baryons, 
respectively. The time evolution of the abundances of 
normal nuclei and X-nuclei are displayed in Figs. [5^ and 
[2b. 

At high temperatures 100 ^ Tg > 10, neutrons and 
protons are the main constituents of baryonic matter 
in the universe since photonuclear reactions dissociate 
bound nuclei at these temperatures. However, as the 
temperature decreases, bound nuclei can form. At Tg ~ 5 



(T 0.4 MeV, t 4 s) X° particles capture nu- 
cleons so that the abundance suddenly decreases 
(Fig. [2]3). The AT" particles first predominantly cap- 
ture protons to form ^Hx- This is because only pro- 
tons and neutrons exist in significant amounts during 
that epoch. Also, neutron capture reactions are hindered 
as explained in Sec. IH C II The ^Hx nuclei then inter- 
act strongly with background neutrons through the pho- 
tonless ^Jix{'n,p)nx transfer reaction so that nx is also 
produced. These nx nuclei then capture protons to form 
^Hx- As heavier X-nuclei are transformed into other 
nuclei with decreased proton number through {n,p) re- 
actions, they lead to the formation of heavier X-nuclei 
via the reaction pathways shown in Fig. [TJ High energy 
nuclear reactions produce abundant ^Hx and slowly in- 
creasing abundances of ^Hex, ''Hex and ^Hex (Fig. [2b). 
The production of ^Hx via the ppx{n,p)'^}Ix reaction is 
also operative. The AT-nuclei increase their nucleon num- 
ber gradually until the temperature decreases to Tg ~ 1 
(T ~ 0.1 MeV, t - 170 s). Then at Tg - 1, a drastic 
increase in the nucleosynthesis of X-nuclei occurs. 

Nuclear reactions at low temperature are important in 
determining the final elemental abundances for normal 
nuclides. Reactions at relatively low temperatures, how- 
ever, are hindered by the Coulomb barriers. Nuclides 
with small atomic numbers, therefore, are more easily 
processed at low temperature. In addition, reactions trig- 
gered by abundant nuclides play important roles in nu- 
cleosynthesis since the reaction rates are proportional to 
the abundances of the reactant nuclei. 

In BBN, the abundance of deuterium is very high at 
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'Bex —/Bex 



ppx ^Hex — ► "^Hex — ^Hex — ► ^Hex 




FIG. 1; Reaction pathways for the X-midei. Sohd arrows indicate the nuclear reactions, while dashed arrows indicate /3 
decays. Arrows are drawn in the direction of positive Q-value. 



Tg ~ 1 (Fig. [2^). Reactions of X-nuclei triggered by 
deuterons are, therefore, efficient at this epoch both be- 
cause deuterium is very abundant and because its charge 
is low. The most important point is that strong pho- 
tonless nuclear reactions to increase mass numbers exist, 
i.e., {d,n) and {d,p) reactions. These reactions help X- 
nuclei to capture more nucleons and become more deeply 
bound. In this way, X-nuclei are processed at Tg ~ 1 and 
heavy X-nuclei up to ^^Cx are produced in individual 
abundances of IQ-^^ < Ax/H < 10"** (Fig.^). 

This model of BBN with the strongly interacting X° 
particles changes an important aspect in SBBN which is 
that the nuclides ^He and ^Li are unstable to the particle 
decay so that they limit the production of heavier nuclei. 
However, if they are stabilized against particle emission 
by binding with an X^, then {d,n) and reactions 
subsequently link '*Hex with ^Ax, ^ Ax with ^Bx, and 
so on. Light element production is thus catalyzed by X^ 
particles. 

Yields of light nuclides from lithium to carbon are sig- 
nificant in this model. The reactions contributing to the 
production and destruction of light nuclides are summa- 
rized in Table IVIIl The last column in the table shows 



the resulting yields for the case of Yx — 10^^ and a 
long lifetime compared to the duration of nucleosynthe- 
sis; Tx » 10'* sec. Nuclides of mass numbers up to 10 
are produced in abundances larger than A/Il= 10^*^. 

The production of nuclides with A > 10 is not sig- 
nificant through the {d,n) and {d,p) reactions. Small 
amounts of ^^Cx and "'^^Cx are produced, however, via 
the (a,n) and reactions. The ^Lix production re- 

action in this paradigm is ^llex{d,n)^Lix while that of 
^Li in SBBN is the ■*He(d,7)^Li reaction. This latter 
cross section is very small due to the associated hindered 
transition rate through an electric quadrupole transition. 
The ""Bex production reaction is ^Lix(p,7)^Bex while 
that of "^Be in SBBN is '*He(3He,7)^Be whose rate is 
smaller than the former reaction because of the smaller 
abundance of the target nuclide and the larger Coulomb 
barrier. The X-nuclides ^Bex, '^'^Bx and ^*Bx are pro- 
duced through {d,n) and {d,p) reactions, while the nu- 
clides ^Be, ^°B and ^*B in the present Galaxy are thought 
to be produced mainly through nuclear spallation pro- 
cesses of heavier CNO isotopes in the Galactic interstel- 
lar medium (49l . 5C , Isil or via the i^-process in supernova 
to produce ^*B 521 l53l| except for non-standard baryon- 
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FIG. 2: Calculated abundances of normal nuclei (a) and X-nuclei (b) as a function of Tg. For this figure we take the X'^ 
abundance to be Yx ~ Nx/nf, = 10~*, and its lifetime is taken to be long tx = oo. We utilize the reaction rates as 
described in the text. The dashed lines in the upper figure (a) correspond to the abundances of normal nuclei in the standard 
BBN model. These are nearly indistinguishable from the solid lines. 
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TABLE VII: Most Important Reactions for the Production 
and Destruction of Light Nuclides 



TiiipliHp 


Tirr^H 1 1 pi" 1 nn 




yield A/H 


•^Liy 


^IIex(d,ri)''Lix 


^LixCd o'l'^Lix 


5 X 10" 


10 


'Bex 


^Lix(p,7)^Bex 


'^Bex(n,p)^Lix 


1 X 10~ 


11 


'Bex 


«Lix(p,n)«Bex 


«Bex(rf,p)^Bex 


8 X 10" 


-9 


'Bex 


*Bex(d,p)^Bex 


^Bex(d,n)i°Bx 


2 X 10" 


-9 


lO-R 
DX 


^Bex(d,n)"'Bx 


i«Bx(d,n)"Cx 


2 X 10" 


-9 


DX 


i«Bx(d,p)"Bx 


"Bx(p,a)*Bex 


1 X 10" 


14 


i2r< 


'>Bex{a,nf'Cx 




3 X 10" 


12 


13^1 


^°Bx(a,pf'Cx 




2 X 10" 
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inhomogeneous BBN |54| . 

An interesting point regarding boron production in 
BBN with X° particles is that the production of ^"Bx 
is preferred over that of ^^Bx- This would imply more 
abundant ^°B than ^^B after decay. There are no 
processes that predict preferential production of ^°B in 
standard processes to synthesize boron isotopes. Galactic 
cosmic ray spallation nucleosynthesis predicts a ratio of 
"B/ioB- 2.5 (e.g., [4i[53,[5H). Supernova nucleosyn- 
thesis in massive stars produces large amounts of ^^B rel- 
ative to "'^^B through neutrino interactions ^^C{v,v'nY'^G 
and '^^C{v,v'pY^Ji in the carbon layers and ^Li(a,7)^^B 
and ^Be(a,7)"'^^C in the helium layers [s^. [53|. 

We note some important differences between BBN with 
strongly interacting particles and BBN with nega- 
tively charged leptonic X^ particles that only interact 
electromagnetically, i.e. no strong interaction. First, 
the formation epoch of X-nuclei in the X^ case begins 
at Tg ~ 5 which is much earlier than in the X~ case 
which begins for Tg < 0.5 [1^. This derives from the 
fact that the binding energies of nuclides and X° parti- 
cles are of the order of ~ 10 MeV which is much larger 
than the electromagnetic binding of nuclides with X~ 
particles (~ 1 MeV). X-nuclei produced in such a high 
energy environment have a better chance to be processed 
so that many species of X-nuclei are produced in consid- 
erable abundance. In the X~ case, on the other hand, the 
temperature of the universe is already so low that when, 
the X-nuclei finally form, they can hardly be processed 
through subsequent nuclear reactions except for the spe- 
cial cases of the resonant ''Bex(p,7)^Bx reaction [55| and 
the X~-catalyzed transfer reaction '^'i{ex{d-,X^)^'Li [s^ . 
As a result, nuclides heavier than the beryllium iso- 
topes are not produced in significant amounts in that 
paradigm. 

This result for the time evolution of the AT-nuclei abun- 
dances was not predicted precisely in the analytic estima- 
tion of previous studies [l^, US] • Dicus and Tcplitz [l9| 
deduced that the stable hadrons would be preferentially 
locked into beryllium. Our result, however, shows that 
many light elements including beryllium are produced 
abundantly and beryllium is not particularly more abun- 



dant than other nuclides. In the study of [l9|, the strong 
photonless reactions {d,n) and {d,p) were not included. 
This is the reason for the large difference between that 
work and our present results for BBN. 

In the present work we have shown that most of the 
X^ particles end up in X-nuclei after BBN in contrast to 
the estimation in f20'|. We note, however, that the X cap- 
ture reaction rates adopted in the present study are only 
approximate and should be calculated more precisely in 
a more thorough quantum mechanical treatment in the 
future. We have shown that strongly interacting X^ par- 
ticles undergo efficient X capture by protons at Tg 5. 
This leads to the subsequent production of ^Hjc, '^Hex, 
^Hejf, and so on. This does not occur through normal 
nuclei like ^H, ■^He and ''He, but through AT-nuclei from 
^Hj!f. This reaction flow was not considered in Al- 
though our result includes some uncertainty in the bind- 
ing energies and nuclear reaction rates of X-nuclei [orig- 
inating from our assumption of the interaction strength] , 
we nevertheless expect that the X^ particles will end up 
in X-nuclei after BBN as long as the interaction strength 
is very large compared to electromagnetic strength. 

The calculated BBN abundances of light nuclides de- 
pend strongly on the abundance, Yx- As an example, 
a series of BBN calculations was carried out varying Yx 
as a parameter with no change in the assumption of large 
Tx- Although the decay hfetime of X° is assumed to be 
much longer than the time scale of BBN > 10^ s, the X^ 
particles are assumed to have been long extinct by now. 
The final abundances of stable nuclides are obtained by 
removing the X'^ from the X-nuclei, Ax, and allowing A 
to decay to a stable nucleus. The interaction of the decay 
products with the remaining A spectator nucleus, and the 
nonthermal nucleosynthesis triggered by high energy de- 
cay products are neglected here. Such effects, however, 
should be studied in the future. 

Figure [3] shows the calculated abundance ratios of 
^Li/H and ''Li/H (solid curves) as a function of Yx- 
The dashed lines correspond to the mean values mea- 
sured in MPHSs, i.e., ^Li/H=(7.1 ± 0.7) x lO'^^ [s^ and 
^Li/H=1.23l^;^^ x IQ-io [H, respectively The ^Li and 
^Li abundances both increase linearly with Yx- However, 
the final ^Lix and ^Lix abundances per A"-particle are 
independent of Yx- As is shown in Fig. [21 the additional 
^Li is produced mainly from ^Lix while the additional 
^Li is from ''Bex. Comparing this with the case of BBN 
with a negatively charged leptonic X~ particle we 
find that both cases prefer the production of ^Li to that 
of ''Li. The difference between the two cases is the rel- 
ative efficiency of X-catalyzed nucleosynthesis. The ef- 
ficiency of ®Li production in the X'^ case is ~ lO'' times 
larger than in the X~ case when the abundance of the 
X particle (i.e. Yx) is taken to be the same. 

In order to check if there is a possibility to solve the 
^Li and ^Li problems in the present AT^-catalyzed BBN 
model, a parameter search was performed over the decay 
lifetime tx and the abundance Yx of the X^ particle. As 
described below, constraints on the (rx, Yx) parameter 
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FIG. 3: Calculated abundances of ®Li/H and ^Li/H as a func- 
tion of the abundance Yx ■ The dashed lines indicate the 
mean values observed in MPHSs. 



space were derived from observational limits on the pri- 
mordial light element abundances. We could not find a 
parameter region which simultaneously resolves the two 
lithium problems. 



B. Constraints on Primordial Light Element 
Abundances 

Deuterium is measured in the absorption spectra of 
narrow line Lyman-a absorption systems in the fore- 
ground of high redshift QSOs. An analysis of well re- 
solved DI Lyman series transitions in the spectrum of 
the QSO Q093H-072 was performed recently [59l|. With 
the newly deduced value of the deuterium abundance, 
they obtained a mean value for the primordial deuterium 
abundance of log(D/H)=-4.55 ± 0.03. We adopt this 
value and a 2 sigma uncertainty, i.e., 

2.45 X 10"^ < D/H < 3.24 x 10"^ (22) 

^He is measured in Galactic HII regions by the 
8.665 GHz hyperfine transition of ''He"'" [60|. A plateau 
with a relatively large dispersion with respect to metallic- 
ity has been found at a level of 3He/H=(1.9±0.6) x 10"-^. 
However, it is not yet understood whether ''He has in- 
creased or decreased through the course of stellar and 
galactic evolution ,'6?|. Nevertheless, it does not seem 
that the cosmic averaged ■'He abundance has decreased 
from that produced in BEN by more than a factor of 2 
due to burning in stars. Because both the observations 
and theoretical predictions of the primordial deuterium 
abundance are consistent with each other, the deuteron 
abundance appears not to have decreased since the epoch 
of BBN. ''He is more resistant to the stellar burning. 
Therefore, its abundance would not have decreased sig- 
nificantly. Thus, we adopt a two sigma upper limit from 



Galactic HH region abundances of 

^He/H < 3.1 X 10"^ (23) 

We do not give a lower limit due to the large uncertainty 
in the Galactic production of ''He. 

For the primordial helium abundance we adopt Y = 
0.2477 ± 0.0029 from an analysis (6^ of the primordial 
mass fraction based upon new atomic physics computa- 
tions of the recombination coefficients for Hel and of the 
coUisional excitations of the HI Balmer lines, together 
with observations and photoionization models of metal- 
poor extragalactic HII regions. We adopt the following 
range for the primordial ^He abundance within the con- 
servative 2 sigma limits of 

0.2419 <Y < 0.2535. (24) 

An upper limit on the ^Li abundance is taken which 
allows for the possible depletion on stellar surfaces of 
up to a factor of ^ 10 above the observed plateau 
abundance of ^U/}i= (7.1 ± 0.7) x lO'^^ This 
upper limit is not larger than the solar abundance of 
^Li/H© = 1.7 X 10""^*^ [64] which represents a typical 
present abundance. We do not take the lower limit on 
the observations as a lower limit for the production of 
^Li by X-particles, because of the current controversy as 
to whether the ^Li observation is actually a primordial 
abundance. Hence, the lower limit of ^Li for our purposes 
is zero. The adopted constraint on the ^Li abundance is 
thus, 

^Li/H < 10-^°. (25) 

An upper limit on the ^Li abundance is taken to be 
6.15 X 10^^" considering a possible depletion of up to 
a factor of ^ 5 down to the observationally determined 
value of ^Li/H= (1.231:^ j^) x 10"!° A lower limit is 
taken from the 2 sigma uncertainty in the same observed 
value. The adopted constraint on the ^Li abundance is 
therefore 

0.59 X 10"^° < ^Li/H < 6.15 x 10"^°. (26) 

We adopt minimum abundances observed in MPHSs 
as constraints on abundances of ^Be, B and C, i.e., 

^Be/H < 10~l^ (27) 
from [6i,[6E[63,[6|[zi, 

B/H<10"i^ (28) 

from [6i,[73] and 

C/H < 10"^ (29) 

from a compilation of observational data by Suda et 
al iZll]. 
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(sec) 

FIG. 4: Contours in the {tx,Yx) plane corresponding to the 
adopted constraints for the primordial abundances. Contours 
for the mass fraction of *He, Y = 0.2419 (red line) and the 
number ratios of ^He/H=3.1 x lO"'^ (green line), D/H=3.24 x 
10"^ and D/H=2.45 x lO"'^ (black lines), ^Li/H^ 10"^" 
(blue solid line) and *^Li/H=(7.1 ± 0.7) x 10"^^ (blue dashed 
lines), ■^Li/H=6.15 x 10~^° (purple line), ^Be/H=10"" (pink 
line), B/H=10-^^ (orange line) and C/H=10"* (gray line) 
are shown. 



C. Observational Constraints on the Long-lived 
Strongly-interacting Particles 

In order to study the effects of X'^ decay we cal- 
culated a series of BBN models in which we varied 
the decay life tx and abundance Yx of the X° parti- 
cle while the baryon-to-photon ratio remained fixed at 
77 = 6.3 x 10^1° [ll. 

Figure 2] shows a contour of the "^He mass fraction 
Y = 0.2419 (red line) in the (txXx) plane. Above 
this contour, Y < 0.2419 which breaks our adopted limit 
on the primordial ^He abundance. Also shown are the 
contours for the upper limits of D/H < 3.24 x 10^^ 
(black line), ^He/H < 3.1 x lO^^ (green line), ^Li/H 
< 10-1° (blue solid line), ^Li/H < 6.15 x 10-^° (pur- 
ple line), ^Be/H < IQ-^^ (pink line), B/H < lO'^^ (or- 
ange line) and C/H < lO"* (gray hue). The contours 
for the lower limit of D/H > 2.45 x 10"*^ are drawn by 
other black lines. The blue dashed lines correspond to 
the observed level of the ^Li abundance in MPHSs, i.e, 
'5Li/H=(7.1 ± 0.7) X 10"^^. The upper right region sur- 
rounded by the contours is excluded based upon the ob- 
servational constrains. 

These constraints in the {tx,Yx) parameter space can 
be summarized as follows. When the decay lifetime of the 
is short {tx < 30 s), no constraint is imposed from the 
observed elemental abundances. When the decay lifetime 
is longer (30 s < rx < 200 — 300 s), the constraint from 
the upper limit on the ^Li abundance is the strongest and 
it implies Yx < 10~* — 10~^. When the decay lifetime is 
longer (200 — 300 s < tx < 2x10^ s), the constraint from 
the upper limit on the B abundance is the strongest and 



it implies Yx < 10~® — 10"^"'^. Finally, when the decay 
lifetime is very long but still much shorter than the age 
of the present universe (2x10'^ s < tx ^ 4 x 10^'' s), 
the upper limit on the ^Be abundance is the strongest 
constraint and it implies Yx < 10"^^ — lO"^^'^. Since the 
relic abundance of the strongly interacting X particles is 
estimated to be Yx ^ 10~® (Eq. [I]), the decay lifetime 
should be tx ^ 200 s. Thus, models suggesting long- 
lived colored particles with tx ^ 200 s are rejected. 

The shapes of the contours on Fig.[5]reflect the forma- 
tion epochs for the X-nuclides (see Fig. [2]) . We are not 
interested in the case of a large abundance of Yx ■ This is 
because the relic abundance of a long-lived SIMP can not 
be very large (see Eq. [1]) although there are significant 
influences of X° particles on light clement abundances 
if the lifetime tx is long. In the lifetime region of 30 s 
< Tx < 200 — 300 s, ^Lix overproduction is the cause of 
the exclusion through the ^Li abundance constraint. One 
can see that the abundance of ^Lix has a peak at Tg '--^ 1 
and is destroyed below this temperature (Fig. ^^Bx 
and ^Bex form at Tg < 1 in almost the same amounts. 
In the longest hfetime case (2x 10^ s < rx 4 x 10^^ s) 
the constraint from ^Be is stronger than that from B, 
since the observed minimum abundance of ^Be is about 
one order of magnitude lower than that of B. 

A new solution to the ^Li or ^Li problems was not 
found in this work. In the reactions included in our BBN 
calculation, there are no reactions which can effectively 
destroy ^Li or ''Be. When a large amount of ^Li is pro- 
duced in this X'^ catalyzed BBN, the abundances of co- 
produced beryllium and boron arc unrealistically high, 
and are therefore excluded. The ^Lix abundance has 
a peak at Tg ~ 1 in Fig. [^b. Even if the X^ particles 
decay during this epoch, the remaining ^Li is easily de- 
stroyed via the ^Li(p,a)'^He reaction which also destroys 
the ^Li produced via the ^He(c?,7)^Li reactions in the 
SBBN. Hence, the preferential production of ^Li is im- 
possible in this model, although corrections for effects of 
X^ decay may resolve this as found in Refs. Q- 



IV. CONCLUSIONS 

We have investigated the effects on BBN of a hypo- 
thetical long-lived strongly interacting massive particle 
(SIMP) The strength of the interaction between the 
X'^ particle and normal nuclei is assumed to be similar 
to that between normal nuclei. Under this assumption, 
binding energies of X particles to nuclei have been esti- 
mated from the Shrodinger equation with a nuclear and 
Coulomb potential. Using these binding energies, the 
reaction Q-values for many nuclear reactions involving 
nuclei bound to X (X-nuclei) were calculated. Reaction 
rates for X capture by nuclei, and the nuclear reaction 
rates of X-nuclei were estimated using information from 
existing reaction rates of normal nuclei. We calculated 
the light element nucleosynthesis simultaneously taking 
into account X capture by nuclides and the nuclear re- 
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actions of X-nuclei along with the standard nuclear re- 
actions. The conclusions are as follows. 

First, some X-nuclides like ^Hex, ^Lix and ^Bex 
are stabilized against particle decay since the binding 
energies of nuclei and X particles are very large ~ 
O (10 MeV) and the total binding energies of X-nuclei 
(i.e. the binding energies of nuclei from separate nucle- 
ons plus the binding energies of X-nuclei from separate 
X particles and nuclei) change significantly from those 
of normal nuclei. Accordingly, there are several reac- 
tions whose Q- values become negative when nuclides are 
bound to X particles. 

Secondly, at Tg ~ 5, the X° particles capture nucleons 
so that the free X^ abundance decreases suddenly. Then 
high energy nuclear reactions produce abundant ^Hx. 
The X-nuclei then increase their nucleon number grad- 
ually until the temperature decreases to Tg 1. When 
the abundance of deuterium becomes high at Tg ~ 1, 
strong photonless nuclear reactions [i.e. {d,n) and {d,p)] 
produce heavier X-nuclei up to ^^Gx- 

Thirdly, from a comparison between BBN with X'^ 
particles and BBN with negatively charged leptonic X^ 
particles, we find that the bound state of X-nuclei forms 
earlier in the case (Tg ~ 5) than in the X~ case 
(Tg ^ 0.5). This is due to the larger binding energies of 
nuclei and X'^ particles. Since X-nuclei are produced in a 
high energy environment, they can be processed to heav- 
ier X-nuclei which are produced in considerable amounts. 

Fourthly, we do not find a solution to the ^Li or ^Li 
problem in this BBN model with the X'^ particles. There 
is a possibility that ^Li is produced in an abundance 
even more than that observed in MPHSs. However, the 
coproduced abundances of heavier beryllium and boron 
nuclides constrain this possibility. 

Fifthly, constraints on the lifetime and abundance of 
the X^ particle are derived. The following constraints 
on the abundance: Yx < 10"^ - 10"^ (for 30 s < 
Tx < 200 - 300 s), Yx < 10"*^ - 10"" (for 200 - 300 s 

< Tx < 2 X 10^ s) and Yx < 10~" - lO^^^-^ (for 2x 10^ s 

< Tx <C 4 X 10^^ s) have been deduced based upon 
the observational constraints on primordial light element 
abundances. These constraints reject models with long- 
lived (tx ^ 200 s) colored particles based upon their relic 
abundance, i.e., Yx ~ 10^*. 



Although we made an assumption about the interac- 
tion strength of the X'^ particle, the magnitude of the 
effect of long-lived SIMPs on BBN would not change by 
more than a few orders of magnitude from the result of 
this study unless the interaction strength is much smaller. 
We also expect that the result would not change by many 
orders even if a SIMP has a charge of ±e. This is be- 
cause the nucleosynthesis of X-nuclei occurs very early 
in the universe when the nuclear reactions are not signif- 
icantly hindered by the Coulomb barriers. More realistic 
estimates of the reaction rates utilizing a quantum me- 
chanical treatment would be necessary to more precisely 
study the effect of X° particles on BBN. Moreover, the 
direct interactions of decay products with the remaining 
nuclei A in the decay of X^ in an X-nucleus Ax should 
be studied in the future in order to better estimate fi- 
nal abundances of the light elements after the X'-' decay. 
Nevertheless, the calculated light element abundances in 
this study should not change by more than an order of 
magnitude from the effects of the decay of X° considering 
the fact that the cross section for the elastic scattering of 
pions by ""^^C is roughly one half of the total cross section 
so that half of the ^■^C nuclei are left intact even when 
effects of the hadronic decay are included [731] . 
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